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We study the possibility of using matter wave interferometry techniques to build a 

O 

gravitational wave detector. We derive the response function and find that it contains 
(3JT), a term proportional to the derivative of the gravitational wave, a point which has 



been disputed recently. We then study in detail the sensitivity that can be reached by 
such a detector and find that, if it is operated near resonance, it can reach potentially 
interesting values in the high frequency regime. The correlation between two or more 
of such devices can further improve the sensitivity for a stochastic signal. 

PACS numbers: 04.80.Nn,04.80.-y,95.55.Ym 

I. INTRODUCTION 

The use of a matter wave interferometer for the detection of gravitational waves has 
been recently proposed in d|. Analogously to what happens with light interferometers, 
matter-wave ones are meant to measure the phase shift induced by the passage of a gravi- 
tational wave (GW); however, dealing with atoms rather than photons implies some crucial 
differences. 
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For instance, what is measured is the effect of the GW on the speed, rather than on the 
position, of freely falling objects (which are the interfering atoms in atomic interferometers 
and the mirrors in the traditional light ones). 

More than that, a matter wave interferometer is sensitive also to the time derivative of 
the GW, h, thus making it convenient to operate it in the high frequency regime. 

The analysis performed in |.L] contains actually some flaws. To begin with, the movement 
of the various parts of the interferometer has not been correctly taken into account for 
frequencies above the mechanical resonance. As a result, the actual sensitivity attains its 
maximum value at the resonance, dropping sharply above it. This means also that the size 
of the device can not be too large, otherwise the lowest mechanical resonance would be too 
low and the sensitivity at that point just not good enough. 

Secondly, as noticed in P , the passage of the gravitational wave affects also the departure 
times of the two interfering paths of the atomic beam, a fact that has not been included 
in the original proposal. Actually the authors of j^J go further and state that, once this 
effect is taken into account, the phase shift does not depend anymore on h, thus completely 
spoiling the interesting high frequency behavior of the response function. We claim that 
this statement is incorrect, because it is based on the crucial assumption that the atoms' 
trajectories of both arms have exactly fixed endpoints, coincident with the detector location. 
In other words, the authors of 2] have analyzed only the case in which the passage of the 
GW does not change at all the extrema of the atoms' trajectories. Such constraint is not 
physically motivated as the typical displacement induced by a GW on the atom position is 
many orders of magnitude smaller than the transverse spread of the atomic wave-packet, not 
to mention the detector acceptance, thus not preventing the atomic beams from interfering 
and being detected. 

We indeed find that, once the fixed endpoint condition is relaxed, a term proportional to 
h is still present in the expression of the phase shift. 

We then explore the potentiality of a GW detector inspired on (an amended version of) 
the proposal of and we require it to be based on present technology. Such a device is 
thus in principle immediately realizable at relatively low cost and in several copies. 

The paper is organized as follows: in the next section we reproduce the relation between 
the dephasing measured at the interferometer and the amplitude of the GW. The discussion 
will be quite detailed, in order to emphasize the points of contact and disagreement with 
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the treatment given in |2(. In section ITTT1 we then examine the principal sources of noise that 
can limit the sensitivity of a typical state of the art interferometer and we plot and discuss 
the relative sensitivity curve. Measurement of static or quasi-static gravitational or inertial 
fields by means of matter interferometers is already a mature discipline: see for instance the 

n n 

pioneering work of [3 J and the recent measurement of the Sagnac effect reported in |4j . The 
extension of this discipline towards the detection of more rapidly varying signals is not a 
trivial task and, as it is a crucial point in our treatment, we will discuss in detail the main 
limitations to a fine time resolution of the signal, and consequently to the detection of high 
frequency GW's. In section HXl we will discuss the possibility of correlating the output of two 
small atomic interferometers and we will compare the obtained sensitivity with the actual 
bounds on stochastic GW signals. Final comments and remarks will conclude the paper. 

In all our equations we will write h explicitly (not to be confused with the gravitational 
perturbation h), while keeping the usual notation for the speed units c = 1. 

II. THE RESPONSE FUNCTION 

Here we compute and discuss the response function for the GW detector based on atomic 
interferometry. 



The spacetime metric in presence of a GW propagating along the h direction can be 
parameterized in the TT gauge as 



with h\ = rfhij = 0. This corresponds to the choice of a free-falling reference frame. 

We find convenient to work in a different gauge, related to the TT one by the coordinate 
transformation 



A. 



General considerations 



.2 



dt 2 + [5ij + hij(t — h ■ x)]dx l dx^ , 



(1) 



t->t' = t 



X = X 




(2) 



which returns the metric 



ds 2 = -dt 2 - dtdx l h ij x j + [5^ + (\i ik x k nj + h jk x k ni)/2]dx l dx j , 



(3) 
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where it is understood that the components of the metric perturbation are the same as 
in the TT gauge. In this coordinate frame the connection coefficients vanish at the origin 
of coordinates, making geodesic equations look like geodesic deviation equations. Moreover 
it can be read from eq. © that the GW does not affect physical distances in the plane 
orthogonal to n. Because of this, in Q] the coordinate frame specified by is named 
geodesic and the (jSJ) rigid frame. 

In the latter, the explicit non-relativistic form (including the next to leading term in the 
velocity) of the action returning the geodesic equation at first order in hij 

V i = -h\X j (l — 2n-V) — -h jk ViV k rj (4) 
2 2 

is then 

S = jj{ v2 - kjV'X^l - h ■ V)) dr , (5) 

with V i = X\ 

Since in a typical interferometer the atom wave function is essentially peaked around the 
two classically allowed trajectories (which will be called the (+) and (— ) paths henceforth), 
the phase in the final point can be calculated, for a given path, as: 

$(/) = *(<) + \s{i -> /) , (6) 

where S(i —>■ /) is the classical action from the initial to the final point. Notice that in 
general one has $(i)( + ) ^ $(i)(~) since the flight times along the two paths are not the same, 
because atoms reaching the detector at the same final time, have entered the interferometer 
at different initial times. 

It is now useful, in order to calculate the action S{i — > f), to introduce the splitting of 
the atoms position and velocity X, V as X = X + 5X, V — Vo + w, where X , Vo are the 
position and velocity in absence of the GW. 

Taking also into account, as pointed out in Q], the shifts Stij of the initial and final times 
tij of the trajectory the action reads as follows 

•tf+stf r _ i . 

-/ 

2 



S ~ S + m 
~ So + m 



u+su 
tf 



V -w--h ij V^Xi(l-h-V , 



Vo-w-^hjVjX^l-n-Vo) 



dr ~ 

(7) 



dT + —V 2 (8tf-5ti 
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where S is the value of the action computed over the unperturbed trajectory. 

Up to first order in the GW perturbation , the initial phase for each path is $(i) 
— ^Yq{U + Sti). This result allows us to rewrite eq. (JHJ) as 



m m 2 J t . 



V .ti- -hyVSXiil -h-V ) 



dr + y 2 (%-^)-^|%(8) 



We now notice that the physical distance between the initial (final) point of the perturbed 
trajectory from the unperturbed one is 

AX U = 8X iif + VoSt ijf , (9) 

and that SXf — 5Xi = f*? wdr, thus giving 

-*(/) = ^ - S + St f ) - \ !\v^Xi{l - n ■ V )dr + V ■ (AXf - AXi) . (10) 
m in 2 2 J t . 

The first two terms do not provide any information about the GW: the first does not depend 
on h, while the second is equal for both paths (since the atoms arrive at the detector 
simultaneously) and thus cancels out when taking the phase difference. We also notice that 
the term proportional to w has disappeared or, better, it has been rewritten in disguise as 
part of the term in AX. 

r |—| 

The analysis performed so far is equivalent to the one done in |2j , whose authors however 

further elaborate on eq. (jl(J|) by setting AX it f = 0. This implies a restriction only to 
such /ijj's that do not change at all both the coordinates of the extrema of the atoms' 
trajectories. However, as stated in the introduction, any position change induced by the 
GW is much smaller than the transverse spread of the atomic wave-packet, and than any 
detector acceptance. This means that the interference and the detection take place also for 
a GW affecting the endpoints of the atoms' trajectories, and the condition imposed in [3] 
does not describe a realistic experimental setup. 

Clearly, by setting to zero the last term in eq. (jlUj) , one loses any possible information on 
h, since the only remaining term, the next to the last, can be calculated by parts and shown 
not to contain any derivative of the GW. This is the ultimate reason of the claim contained 
in about the absence of a term in h in the expression of the phase shift. 

It is useful to end this subsection by noticing that one component of AXij can be 
expressed in terms of the other by means of the direct substitution of 5tij in eq. (jOJ), thus 
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(a) 



(b) 



Figure 1: Schematic configuration of the interferometer. Choice of the polar angles is also displayed. 



giving 



AY f - AYi 



w y dr 



V v 

yx 
1 



AX f - AX 



w x dr 



(11) 



B. The Mach-Zender configuration 

The exact form of the response function depends on the details of the experimental 
setup. For the sake of concreteness, we will consider the Mach-Zender configuration which 
is displayed in Fig^ which we assume to be realized by means of properly chosen diffraction 
gratings. Any other configuration and/or other deflecting device (mirrors, lasers) can be 
studied along the lines which we are going to expose. In the absence of GW, the classical 
upper and lower trajectories are 



x (±)(i~) = V rca, 



F (±) (r) = ±V rsa for < r < T/2 

y ( ±)(r) = ±V (T - r)sa for T/2 < r < T . 



(12) 



with ca = cos a, sa = sin a and we have from now on tj = 0, tf = T. We can now explicitly 
evaluate eq. (|1U|) along each arm of the interferometer, that is from O to A, from A to C, 
do the same for the (— ) path and then take the difference. With the help of eq. the 
following expression is obtained: 

"T r n/y\2 i 



-A ± $(C) = A : 
m 



v, 



o 



-w x + V*w v - -MM(1 -n-V ) 



dr+ 



+^A ± (AX C - AX ) 



(13) 
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where A^ means that the difference between the upper and the lower path is taken. The 
first and the second terms in the integral and the last term are the ones discarded by Roura 
et al, and in particular the parts containing w x,y , which come from the expression of AY 
through eq. (JUJ, will give a contribution proportional to h. 

As to the last term, we observe that A ± (AXc > o) = 0, that is the variation of the x- 
coordinate of the extrema of the trajectory are the same for the upper and the lower path 
and we can then discard the last term in eq. (fT3|) . 

From now on we restrict to the case of normal incidence, n« = 5i Z ; in this case the symmet- 
ric geometry shown in Fig^ makes the device sensible only to the polarization represented 
by h xy = h x , because the h + polarization affects in the same way the (+) and the (— ) paths. 

The dependence of w x,y on h x can be derived by expanding to first order in h the geodesic 
equation (jlj) around the unperturbed trajectories (fT2^) . 

In doing this one has also to take account of the junction conditions at the grating. We 
model their effect as a rotation the atomic velocity by an angle 2a Q|: 

V {±) (T+/2) - V a ,b = 1Z[T2a] ■ {v {±) {T~ /2) - V a ,b) , (14) 

where we have denoted by /2 the instant of time immediately after (before) the in- 
teraction between the grating and the atoms. The formula accounts for the velocities Va,b 
acquired by the grating placed in A and B because of the passage of the GW; these terms 
also contain the derivative of the GW and have been inaccurately taken into account in 
Q]. Indeed we will see that their contribution suppresses the sensitivity in the very high 
frequency regime. 

By expanding to first order in h one obtains the junction conditions for w 
w* (±) (T+/2) = c(2a) wl ±) {T-/2) + (1 - c(2a))V^ ± (wf ±) (T-/2) - V\ B ) s{2a) , 

(15) 

w y {±) (T+/2) = c(2a) wf ±) (T-/2) + (1 - c{2a))V% B =F («f ±) (T-/2) - V^) s(2a) . 
After some calculations, one arrives at the following expression for the dephasing: 



^-A ± $(C) = V satga{2\ sn ( /-/ )/,.,.{ -jdr-Kf 




T 



;ie) 



The square brackets in (|T^|) contain terms involving, both explicitly and implicitly, the 
time derivative of the GW; this is the most important contribution at frequencies / such 
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that T -1 <C f, and it is all due to the interaction with the central grating, which introduces 
a discontinuity in w as an effect of the junction conditions (fT3j) . In particular, the term 
involving V% B represents the direct effect of the state of motion of the grating, which imparts 
a kick to the atoms. We will show shortly that this effect cancels the preceeding h x {T/2) 
contribution in the limit f ^> fo, being fo the typical resonant frequency of the apparatus. 
Such frequency is roughly given by f = c s / (2L), where c s is the sound speed of the material 
the experimental apparatus is made of, and L its linear size (of order VqT/2 in our case). 

In order to obtain a closed form for the phase we need to express the velocities of the 
central grating with respect to the point O in terms of known quantities. This can be done 
by solving their geodesic equations, which are damped harmonic oscillator equations with a 
source term given by the GW: 

*A,B + —^Xa,B + ( 27T fo) 2 ^A,B = 2^) X A,B , ( 17 ) 

where Q is the quality factor of the oscillator and Xa,b the position of the grating. Being 

h(t) = J dfh(f)e- i27rft , h(f) = J dt h(t)e l27Tft , (18) 

our convention for the Fourier transform, the perturbed velocities of the gratings can be 
easily written as 

vj - v. s - MS _ xi) = , (19) 

where in the Fourier representation of the solution we neglected the transient term, i. e. the 
motion of the apparatus not induced by h, which is equivalent to assume that it was at rest 
before the arrival of the GW. 

We can finally write the result for the Fourier transform of the dephasing. Here we report 
only the terms which dominate in the regime T" 1 /, where the matter- wave interferometer 
is more performant: 

A«* (/ ) = V1 A tgV (, _ 1 _ /oV/2 1 +;/o/(/Q) ) S, , (20) 

where A = (VqT) 2 sa ca/2 is the area of the interferometer. 

We now see explicitly that for f/fo ^> 1 the h x effect vanishes. This was to be expected 
as in this limit both the atoms and the diffraction gratings are freely falling and during 
the deflection of the atoms by the experimental apparatus, assumed instantaneous, the 
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GW has no net effect. Typical time scale of such interactions are smaller than the usee, 
so indeed instantaneous in our problem, see Sect lHIl We underline that by using other 
configurations and/or deflecting devices, a response function quantitatively different from 
(120 j) can be obtained. For instance assuming that the atoms are deflected by bouncing off 
mirrors (which turn out to be the most efficient devices, see sec JIII C\ instead of by grating 
rotation at A and B, the prefactor of h would change from Atg 2 a to A, and dependence 
on h(0) — h(T/2) might be introduced; however the response function will have the same 
qualitative properties, that is the presence of a term proportional to fh , which is however 
suppressed for / ^> f . 



C. Generic incidence direction 

For a GW impinging on the interferometer from a generic direction 

Q = (sin 9 cos 0, sin 9 sin 0, cos 9) (21) 

eq. ()20|) needs to be generalized. The two possible polarizations can be written in terms of 
two mutually orthogonal unit vectors m and n as 

e ab = ™arhb ~ n a n b , 
e ab = mam + n a rh b . 

The unit vectors m, n must lie on the plane orthogonal to Cl and must be mutually orthog- 
onal; one possible choice is 

m = (cos 9 cos (f>, cos 9 sin 0, — sin 9) , ^ 
n = (— sin 0, cos 0, 0) . 

Keeping the axes as in fig^the interferometer will still be sensitive to the h xy component 
of the GW, but for generic incident direction Cl one has h xy = F + h + + F x h x , where the 
pattern functions F +jX are given by 

1 + cos 2 9 . n , 

F+ = sin 20, 

2 (24) 

F x = cos 20 cos 9. 

With respect to a different pair of axis m', n', rotated with respect to m, n by an angle 
ip in the m,n plane (it is enough to have < 0> < vr/4), the F + x of (|2*4*j) becomes F' + = 
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cos2^F + + sm2ipF x , F' x = cos2x(jF x — sm2ipF + . The simple case reported in eq. (|2*0|) 
corresponds to the choice 9 = (f> = ip = 0. 

Strictly speaking this is not the only effect which has to be taken into account to describe an 
oblique incoming direction, as the action (JJJ) involves a term proportional to h ■ V Q , which is 
not vanishing in the case of not normal incidence. Anyway this effect represents a relativistic 
correction to the leading term and thus it is negligible. 



III. ANALYSIS OF THE NOISE 



Let us now concentrate again on the case of ideal incidence angle and polarization. The 
sensitivity of this kind of interferometer is determined by all the sources of noise that enter 
the actual measurement of 



A. Shot noise 



The possibility to detect GW's with an atomic interferometer is fundamentally limited 
by phase/particle-number uncertainty, and so by readout shot noise in the atomic detector. 
In the case of atomic interferometry it is a fair assumption that the total number of atoms 
detected per unit of time is perfectly known. In that case 0, assuming contrast and detection 
efficiency equal to 1 and a low density beam, the error on phase measurement is inversely 
proportional to the square root of the number of detected atoms N 

a| hot = -U = , (25) 

where T is the atomic flux and % the averaging time. If the phase shift to be measured 
is time independent, in principle % can be increased (collecting more atoms) to reduce the 
error on the phase. However if the phase shift is time dependent, it is useful to perform 
a discrete Fourier analysis. To detect a variable $(£) we have to bin the data over a time 
interval % = At thus obtaining a row of data $(tj). We assume that the experiment is 
calibrated so that the number Nj of atoms detected at each bin is given by 

Nj = FAt(l + • (26) 
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Then each of the is a random variable, with average zero (in the absence of signal) 



and variance <jj — 1/ y\FAt. The Fourier transform of the data reads 

M 

$(/*) = At ^e- i2 ^ M $(t,), (27) 
i=i 

where the discrete frequencies are fk = k/t ex , with k in the range 1 < k < M, being M the 
number of bins, and we also defined the total experiment time t ex = MAt. Using the central 
limit theorem we can state that the $(/fc)'s are Gaussian distributed random variables with 
vanishing mean and variance given by the sum in quadrature, thus obtaining (in the limit 
M> 1) 



for both the real and the imaginary part of We note that grows with the integrated 

time of the experiment (we are summing, rather then averaging, when we take the Fourier 
transform) and it is frequency independent. The latter property may seem counterintuitive 
but one has to remember that even if higher frequencies are characterized by higher relative 
errors on a range of time tk = l/fk, nevertheless fits more times in the total experiment 
time and the two things compensate. Notice that in the previous analysis of the shot noise, it 
was implicitly assumed that the binned phase was constant over the binning time. However 
the result is unchanged even if the previous assumption is relaxed j^. 

If we now plug eq. (J2HJ) into the relation (ffllj) . to be more precise, between h(f) and $(/), 
we can estimate the minimum detectable signal as 

h ■ (f) V_WZ (291 

where Fh(f) can be derived from eq. (|T6|) 

pW /Tp/ n _, 1 ( sin (tt/T/2) \ 2 

hU) ~ l-Wf 2 + ifo/(fQ) V «fT/2 ) ■ (30) 

We note that Fh goes to zero in both the low and high frequency limits, f'T — ► and / > /o, 

it is of order 1 for 1/T < / < / and attains its maximum at resonance, where Fh ~ Q. 

To make contact with common literature on GW detectors it is convenient to express the 

sensitivity in terms of h c (f) = h m i n (f)/y/%^. h c has dimension of square root of time and 

it is measured in Hz _1//2 ; the main advantage of its definition is that it does not depend on 

the observation time t Px . 
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B. Other noise sources 



As previously stated, the matter-wave interferometer reaches its best sensitivity at / ~ 
fo = c s /{VqT) = c s /(2L), while for f > fo the sensitivity is drastically reduced because in 
this regime the central grating is not anymore rigidly attached to the beam splitter and is 
thus practically comoving with the atoms. 

Actually, the maximum detectable signal frequency f max is related to the binning time 
At through / max ~ 1/ At. Then, for a given At, we must ensure that the flux is large enough 
in order to have enough atoms to arrive in each bin, that is we must require T ^> fo- 

Moreover, we remind that in the previous section we neglected the spatial dependence 
modulation of the GW. This can be done in the long wavelength approximation, which 
ceases to be accurate when the wavelength of the GW is of the same order of magnitude 
of the length of the atomic paths. This enforces the condition gives JqL <C 1, being L the 
total length of the interferometer; this condition is always satisfied because of the smallness 
of the sound speed c s with respect to the speed of light. 

Even with these qualifications, beside the shot noise there are other possible noise sources 
that should be taken into account. 

First of all, since the atomic beam is prepared with a finite velocity spread Svq, then there 
is the possibility atoms entering the apparatus at the same time arrive at the detector at 
different time bins just because of this spread, thus introducing spurious noise in the signal. 
This leads to a minimum binning time t m j n ~ T5Vo/Vq which thus enforces the condition 
(here and henceforth we systematically neglect factors of order 1) 



that is the influence of the thermal vibrations of the gratings (or mirrors) on the atoms' 
movement. 

By a straightforward use of the the fluctuation dissipation theorem we can compute the 
(single sided) power spectrum of the velocity of the deflector S v as 



being Z(cu) the impedance function, i.e. the function relating the Fourier transform of the 
force and the velocity according to F(lu) = Z(uo)v(uo). The explicit form of the impedance 





(32) 
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can be deduced from (JT7j) and (fT§|) and it gives 

S = 4kBT U2UJ ° /Q (33) 

V M dcf (0,2 _ „2) 2 + (^/Q) 2 ' 

being Md e f the deflector mass. 

This induces a noise power spectrum in the phase, which can be estimated according to 
eq. (fTBj) as 

S<s>, v = (-tVqT saiga*) S v . (34) 



To keep under control the thermal noise we must require that S$ tV < Sq, :S h ot = 1/2.F; given 
the shape of the spectra, it is enough to impose the constraint at resonance, thus obtaining 
the following estimate: 

T < 10K ( Md A (Q-} 1 ( - \ 1 ( Wo ) ( L \ 2 (— \ 2 f35) 
VlKgyVlO 5 / Vl° 8 s~V V4- 10 5 Hz/ V0.15my \Am p ) ' 1 ' 

which has been computed for a = tt/A. This means that the apparatus can not work at 
room temperature; on the other hand the required temperature is not very severe and can 
be reached with standard cryogenic techniques. 

Other kinds of thermal noise do not give more severe conditions than (|35j) . The atomic- 
thermal shot noise (deformation of the deflector because of the heating caused by the im- 
pinging beam) is completely irrelevant here because of the low power of the atomic beam. 

The Brownian motion of the deflector surfaces (structural damping) and the thermoelastic 
damping can be estimated along the line of The resulting induced spectral density of 
the position of the deflector surface S x has to be compared with S v /u 2 derivable from (J3*3~j) . 
With different numerical values for the thermo-mechanical specifications of the deflector 
(Young modulus, mass density, specific heat, expansion and diffusion coefficients) and an 
atomic beam size of 1mm we find that at resonance and at a temperature of 10K these noises 
are several orders of magnitude smaller than the noise induced by S v . The noise induced 
by the structural damping is competitive with the noise induced by Brownian motion of 
the deflectors characterized by eq. ()33)1 only well outside the resonance (at a frequency three 
order of magnitude smaller than /o), where the interferometer sensitivity is anyway degraded 
by the peculiar behavior of the response function (|29jl. 

Finally there are other effects, related specifically to the details of the interferometric 
apparatus used, which determine a reduction of the contrast of the interferometric figure 
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and of the atomic detection efficiency. They correspond to the replacement of eq. (|26|) with 
the expression Nj = 7?JFAt(l + 5g^) where rj is the detection efficiency and C is the contrast. 
Consequently, we can take these two effects into account by multiplying eq. (J2*nj) by a factor 

v/c 

C. Technical specifications and best sensitivity 

We now compute a realistic estimate of the sensitivity of an atom interferometer. Its 
size is a crucial factor, as the sensitivity increases with the area but at the same time it 
increases with the signal frequency and the maximal frequency /o is inversely proportional 
to the length of the arms L. 

With L = 15 cm and using the sound speed of graphite (c s = 18350m/s) we can have /o — 60 
kHz. It is also important that the opening angle of the interferometer is not too small; such 
a parameter depends on the technology employed. 

If nanofabricated diffraction gratings are used one has tana ~ 2h/mVoa with a the grid 
spacing. Typical values used in interferometers are a ~ 200nm giving tana ~ 10~ 4 for 
heavy atoms (Vo ~ 10 3 m/s, m ~ 200m p , being m p the proton mass). However gratings 
with smaller spacing up to a = 32nm have been built Q, as pointed out in [jj. 

Reaching a large a with optical beam splitters is instead challenging, since one has to rely 
either on very small wavelengths or on high order Bragg scatterings. In the latter possibility, 
very well collimated beams are required and in general intensity is strongly reduced. 

Using mirrors large angle deflection can be achieved allowing to build a squared interfer- 
ometer. At present such mirrors have been built using crystal surfaces |s| and utilized with 
Helium atoms; it follows that for our device it is advisable to use Helium atoms, for which 
fluxes T ~ 10 8 sec _1 can be reached. The choice of light atoms like Helium is also necessary 
to overcome crystal surface imperfections ,9]. Given the current ability on preparing crystal 
mirrors, it is unlikely that heavy atoms like Cesium will be used with crystal mirrors in 
a near future Q|. For helium we can take a value for the atoms speed of 10 3 m/s, with 
SVo/Vo ~ 10~ 3 , which ensure the validity of (J3TJ). 

Finally, the sensitivity at the resonance can be improved by choosing a high Q for the 
oscillator: this means paying a price in terms of bandwidth of the detector, but we will see 
that an high quality factor is necessary for obtaining an acceptable sensitivity. Actually 
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matter interferometers are intrinsically wide band, but they can be turned into narrow band 
detectors in order to increase their peak sensitivity. 

By building an interferometer with the previous specifications, one can make it work 
efficiently at the resonance frequency /o — 60 kHz, where it attains its peak sensitivity: 



™. ^ ^li^ ~ 2 • lO-^Hz" 1 / 2 x 



(cm\ ( t y 1/2 ( l y l f^y 1 f Cs y 1 

V 10 J V10 8 s-V \4:m p J V0-15my \W 5 J \1.8 • 10 4 m/s/ 
for a = 7r/4 and m p the proton mass. 

This sensitivity can be compared with the one of a light interferometer at the same 
frequency (at the frequencies we are interested in, h^ lGO is dominated by the shot noise). In 
Fig|2]we compare the LIGO sensitivity curve, as it has been recently reached [f| in the S3 
run, with the sensitivity curve for a matter-wave interferometer endowed with L = 15cm, 
T = 10 8 atoms/s, tana = 1, and a total efficiency of 10%. An extrapolation to high 
frequencies of the LIGOl sensitivity goal is also showed. 



IV. STOCHASTIC SIGNAL: CORRELATION OF TWO INTERFEROMETERS 

For a given h c one needs a much different gravitational wave amplitude according to 
whether the signal to be detected is a burst or a monochromatic wave. Let us consider e.g. 
a burst of amplitude hp and duration time tb, its Fourier transform is hs — h^T^ thus 
giving h c ~ liBy/TB- So, in order to be detectable by an interferometer with sensitivity h c a 
burst must have an amplitude such that 

h B >^- (37) 



This has to be compared to a periodic signal, say a cosine of frequency fp and amplitude 
hp. The Fourier transform in this case gives h(fp) ~ hpt cx , thus the signal is detectable if 

hp>^=, (38) 



and since t ex ^> Tp, one has that periodic signals are easier to detect than bursts. 

Actually the dependence of the minimum attainable h, for a given h c , on the nature of 
the signal is even stronger if correlation of detectors is taken into account. 
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Figure 2: Thick solid lines: sensitivity curve for a Helium interferometer with L = 0.15m, a = 45°, 
fo = 60kHz, Q = 10 5 , J- = 10 8 atoms/s, total efficiency of 10%, and sensitivity curve for the S3 
run of LIGO 6]. Solid line: sensitivity goal for LIGOl extrapolated to the MHz. Dashed line: 
sensitivity curve for two correlated matter interferometers. Dash-dotted line: BBN bound on a 
stochastic background of GW. 

In the case of a GW of stochastic origin a signal even lower than the noise level could be 
detected, provided one has two detectors and their outputs are correlated. We report here 
only the relevant formulas, addressing to [10] for a detailed and systematic treatment. By 
analyzing the output of two detectors for an observation time t ex , and assuming that both 
the stochastic signal and the detector sensitivity are slowly varying on a frequency range Af, 
one can reduce the minimum detectable amplitude according to the following expression 



2,1 



(39) 



which, for t cx = lyear and Af ~ 1Hz, gives an amplification factor of about 10 2 . Corre- 
lation of two light interferometers cannot be realized at high frequencies (/ > 50Hz), as 
the maximum frequency of a signal which can be possibly correlated is given by the inverse 
distance of the two detectors, /o = 6 • 10 4 Hz corresponding to a separation of roughly lKm. 
Of course correlation of a light interferometers with other kinds of devices is possible. 
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Given the moderate cost of a matter-wave interferometer, the construction of several 
matter-wave interferometers is an affordable task: the correlation of n such detectors would 
increase the sensitivity to a stochastic signal by a factor ~ n 1 / 2 ll] with respect to the 
n = 2 case. 

In the 50 kHz region the only known candidate for such a signal is a background of 
cosmological origin, whose strength is limited by the so-called nucleosynthesis bound. The 
amount of gravitational waves produced in the early Universe is constrained by the fact 
that it can not carry too much energy density, as otherwise it will modify the production of 
elements during nucleosynthesis, which is predicted by the standard theory with a remarkable 
accordance with observations. 

This gives a bound on the fraction on energy density that can be carried by these relic 
cosmological gravitons: 



where h is a number of order one that parameterizes our ignorance about the exact value 
of the Hubble parameter. Since the energy density fraction in GW's and h c are related by 



one has that such a cosmological background is out of the reach even of the combined 
operation of two small matter wave interferometers. As this bound on h c goes as f~ 3 ^ 2 
at high frequency, it is useless to push the matter-wave interferometer at extremely high 
frequency in the search for a background of gravitational origin. 

There is a qualification to be done on the last statements: in deriving eq. (}4*Tj) it is 
implicitly assumed that the spectrum of the stochastic signal is smooth enough; the bound 
can thus be avoided if a narrow peak is present in the spectrum. 

V. CONCLUSIONS AND REMARKS 

We have studied the feasibility and the performance of a table-size matter-wave interfer- 
ometer for the detection of gravitational waves. 

We find that its sensitivity in the high frequency regime is proportional to the time 
derivative of the gravitational wave, contrary to what is claimed in [2], but still this is 



6 



(40) 




(41) 
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not enough to compete with the large light interferometers in their frequency bandwidth. 
However matter-wave interferometers can complement their light counterpart in the high 
frequency region (above 10 kHz). 

We have thus studied in detail the sensitivity of a matter wave interferometer around its 
maximal working frequency, and we have found that the atomic beam quality and opening 
angle are the key factors for its sensitivity. Increasing the size improves its sensitivity by 
an area factor, but it decreases the frequency of highest sensitivity. Peak sensitivity can 
be improved exploiting a resonance mechanism, thus making matter interferometers narrow 
band, but it is not excluded to use them as wide band detectors once technological advances 
(atom lasers for instance) will allow to lower noises. 

The correlation of two or more matter-wave interferometers can improve the sensitivity 
for a stochastic signal, but the nucleosynthesis bound makes the detection of a cosmological 
background still out of reach. 

In the frequency region above 10 kHz the only astrophysical sources could be coalescing 
primordial black holes, whose mass can be smaller than the Chandrasekar limit by orders of 
magnitude. However the strength of the gravitational radiation they will eventually emit is 
too small in this frequency range and the signal is decreasing with the wave frequency and 
increasing with the mass of the coalescing objects, thus being stronger in the LIGO band 
rather than in the matter interferometer one. 
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